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On the tidal interaction of a solar— type star with an orbiting companion: 
Excitation of g mode oscillation and orbital evolution 
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ABSTRACT 

We calculate the dynamical tides raised on a non-rotating solar-type star by a 
close stellar or planetary companion. Dissipation arising from a turbulent viscosity 
operating in the convection zone and radiative damping in the radiative core are 
considered. 

We compute the torque exerted on the star by a companion in circular orbit, and 
determine the potentially observable magnitude of the tidally induced velocity at the 
stellar photosphere. These calculations are compared with the results obtained by 
assuming that a very small frequency limit can be taken in order to calculate the tidal 
response (equilibrium tide). For a standard solar model, the latter is found to give a 
relatively poor approximation at the periods of interest of several days, even when the 
system is far from resonance with a normal mode. This behavior is due to the small 
value of the Brunt- Vaisala frequency in the interior regions of the convection zone. 

It is shown that although the companion may go through a succession of resonances 
as it spirals in under the action of the tides, for a fixed spectrum of normal modes its 
migration is controlled essentially by the non-resonant interaction. 

We find that the turbulent viscosity that is required to provide the observed 
circularization rates of main sequence solar-type binaries is about fifty times larger 
than that simply estimated from mixing length theory for non-rotating stars. We 
discuss the means by which this enhanced viscosity might be realized. 

These calculations are applied to 51 Pegasi. We show that the perturbed velocity 
induced by the tides at the stellar surface is too small to be observed. This result is 
insensitive to the magnitude of the turbulent viscosity assumed and is not affected 
by the possibility of resonance. For this system, the stellar rotation and the orbital 
motion are expected to be synchronized if the mass of the companion is as much as 
one tenth of a solar mass. 
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1. Introduction 

Theoretical analyses of the tidal interaction between close binaries can be classified according 
to whether an equilibrium tide is assumed or the dynamical tide is taken into account. The theory 
of the equilibrium tide is based on the assumption that a star subject to the tidal disturbance of 
a companion instantly adjusts to hydrostatic equilibrium (Darwin 1879). A calculation including 
the dynamical tide takes into account the fact that gravity or g modes can be excited in the 
convectively stable layers of the star and that resonances between the tidal disturbance and the 
normal modes of the star can occur (Cowling 1941). So far, dynamical tides have been studied 
only in massive close binaries, which have a convective core and a radiative envelope (Zahn 1975, 
1977; Savonije & Papaloizou 1983, 1984, 1997; Papaloizou & Savonije 1985, 1997; Savonije, 
Papaloizou & Alberts 1995). 

In this paper, we examine the effect of dynamical tides excited by a companion on a solar-type 
star, in which a radiative core is surrounded by a convective envelope. 

This is of particular interest in connection with circularization of solar-type binaries. It has 
been proposed that circularization occurs through the action of turbulent viscosity, originating 
in the convective envelope, on the tide. However, according to Claret & Cunha (1997) (see 
also Goodman & Oh 1997), who have used the equilibrium tide formalism of Zahn (1989), the 
circularization rate resulting from this mechanism is too small by about two orders of magnitude 
to account for the circularization timescales required on the main sequence. 

The tidal response calculation undertaken here is also of interest in connection with the newly 
discovered planets, some of which are found to orbit around solar-type stars with a period 
comparable to that of the high order g modes of the star. One such example is 51 Pegasi (Mayor & 
Queloz 1995; Marcy k Butler 1995). 

In these binaries, g mode oscillations are excited by the companion in the radiative region beneath 
the convective envelope. They become evanescent in the convection zone where they are damped 
by their interaction with the convective eddies. This dissipation leads to an exchange of angular 
momentum between the star and the orbit if the stellar rotation and the orbital motion are not 
synchronized. Here we assume that the orbital frequency is initially larger than the rotational 
frequency of the star. Then tidal interaction results in the decay of the orbit and the spin up 
of the star. If the mass of the secondary companion is considerably smaller than that of the 
primary, the timescale for orbital decay is smaller than the stellar spin up timescale, and the 
companion eventually plunges into the primary. But if the mass of the companion is large enough, 
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synchronization may occur before the binary has merged, stopping further orbital decay. Estimates 
based on the theory of the equilibrium tide (Rasio et al. 1996; Marcy et al. 1997) suggest that 
the orbital decay timescale and the stellar spin up timescale for a system like 51 Pegasi are longer 
than the inferred age of the primary if the companion is a Jovian like planet. 

In this paper we examine the effect of resonances on these timescales, and determine the 
potentially observable magnitude of dynamical tides at the photosphere of a solar-type star. We 
also compare the dynamical tide calculations with the results of an asymptotic analysis we carry 
out in the limit of small frequencies which should correspond to the adiabatic equilibrium tide 
theory. 

The paper is organized as follows: In we study the tidal response of the star to the 
perturbation by a companion in a circular orbit with a period in the range 4-13 days. In § 2.1, we 
first consider the linear adiabatic response and then, away from resonance with a g mode, extend 
the analysis using first order perturbation theory to calculate the torque due to dissipation in the 
convective envelope. This mechanism is then found to be more important than non-adiabaticity 
arising from heat transport in the radiative interior (i.e. radiative damping). However, this is 
not the case in the vicinity of a g mode resonance. There we also calculate the torque due to 
non-adiabaticity in the radiative core using a WKB treatment of the non-adiabatic terms. We 
find that the torque at effective resonances is mainly determined by radiative damping. An 
analysis valid for very low frequencies [equilibrium tide) is given in § |2.2| . 

The orbital circularization timescale for systems initially in non-circular orbits can be derived 



from the response calculations for companions in circular orbits. This is done in § 2.3. We then 



discuss how this might be used to calibrate the magnitude of the turbulent viscosity required to 



fit the observations in S 2.4 



Numerical calculations are presented in §|3| The results assuming the equilibrium tide are given 



m 



34]. In § we present the results of the dynamical tide calculations. We describe the tidal 
response of the star to a companion in circular orbit, give the induced velocity at its surface and 
the tidal torque. We describe the resonances and show that, for the periods of interest of several 



days, they are not expected to affect the orbital evolution of the binary. In § 3.3 we compare 
the calculations based on the dynamical and equilibrium tides. We find that, for the standard 
solar model, at the orbital periods of interest, because of the long timescale associated with 
convection, the equilibrium tide calculations give a relatively poor approximation to the results of 
the dynamical tide calculations. 

We find that the viscosity that is required to provide the observed circularization rates is about 
50 times larger than that simply estimated from mixing length theory and discuss the means 
by which this viscosity might be enhanced in § |3.4| . However, we note that the strength of the 
resonances for orbital periods larger than ~ 8 days and the perturbed velocity at the surface of 
the star are insensitive to the magnitude of the turbulent viscosity assumed. Only for periods ~ 4 
days is the strength of the resonances decreased by a factor ~ 4. The observable width of the 
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resonances as well is reduced when the viscosity is increased. We also give the relation between 
the orbital evolution, circularization and spin up timescales and the orbital frequency in § |3i 



Finally in §0 we discuss and summarize our results, applying them to 51 Pegasi in §4.1 



2. Tidal response to a companion in circular orbit 

The calculations presented in this section will be applied to close binary systems where the 
primary is a solar-type star and the secondary a stellar or planetary companion. The orbital 
periods of interest lie in the range 4-13 days. The rotational angular velocity of the primary is 
assumed to be small compared to the orbital frequency, so that it can be neglected. Quadrupolar 
tidal forcing thus occurs through potential perturbations with periods in the range 2-6.5 days. 

When calculating the tidal response well away from a condition of resonance with a g mode, 
we firstly calculate the tidal response assuming it to be adiabatic throughout the star. First-order 
perturbation theory is then used to calculate the dissipation occurring in the convective envelope. 
The idea here (as is borne out by the numerical results) is that although short wavelength 
g modes are excited in the radiative core, when they are away from resonance they do not play an 
important role in comparison to the global component of the tidal response. Also the variations 
in the convective envelope occur on a comparatively long length scale, making the adiabatic 
approximation a reasonable one. 

When there is a resonance with a high order g mode, the response becomes one with a 
very short length scale such that non-adiabaticity in the radiative core cannot be neglected. 
However, the modes are of high order such that a WKB treatment of the non-adiabatic effects 
is possible and this is used close to resonance where the normal mode dominates the response. 
Such non-adiabatic effects turn out to be more important than the action of turbulent viscosity 
in the convective envelope, with the torque at significant resonances being determined mainly by 
non-adiabatic effects. 



2.1. Linearized equations 

2.1.1. Adiabatic perturbations 

The linearized momentum, mass, and energy equations governing the adiabatic response of 
the non-rotating star to the perturbing potential may be written (see, for example, Unno et 
al. 1989) 



dt 2 



1 d 

-VP' + -7T VP - V*Q 



(1) 
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d = -V.(/«), (2) 
AS = 0, (3) 

where P is the pressure, p is the density, S is the entropy, £ is the Lagrangian displacement vector, 
A denotes the Lagrangian perturbation and the primed quantities are Eulerian perturbations. We 
use the Cowling (1941) approximation, applicable to stars with high central condensation, which 
neglects the perturbation to the stellar gravitational potential. We also have the thermodynamic 
relation 

where T is the temperature and Ti and T3 are the adiabatic exponents of Chandrasekhar. This 
relation together with equation (||) leads to 

d P' 

(5) 



where 



. dlnp 1 dlnP N 2 

A =—, r T~ = ' v 6 ) 

dr 1 1 dr g 

with sgn (N 2 ) x ^/|iV 2 | being the Brunt- Vaisala frequency and g the acceleration due to gravity. 

Following Cowling (1941), only the dominant quadrupole term is considered in the perturbing 
potential due to the companion. For a binary system with a circular orbit, this is given in spherical 
polar coordinates (r, 9, <p) by the real part of 

* T (r, 9, if, t) = fr 2 Y n>m (9, <p) e~ im ^\ (7) 

where the spherical harmonic 



y n , m (9, if) = PH (cos 0) e im * 

with n = m = 2, P^ being the associated Legendre polynomial with indices n and m. Here uj is 
the orbital angular velocity, / = —GM p /4D s , where D is the orbital separation, and M p is the 
mass of the companion. We adopt the same angular and time dependence for the perturbations, 
so that P', d and S' are proportional to Y n>m (9, if) exp (— imuit). The corresponding expression 
for the Lagrangian displacement is 
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Yn^^y ■ (8) 

The factor Y nm {0, <p) exp (— imut) will be henceforth taken as read, so that hereafter the 
perturbations will be taken to depend only on r. Physical perturbations are then found by taking 
real parts after inserting this factor. 

The horizontal displacement ^ is given by the non-radial equation of motion (|l|) : 



06 sin Ooip 



th = ^r{- + fr 2 ). (9) 

m z uj z r \ p J 

This relation together with equation (||) allow P' and p' to be eliminated from the system (Q)-(^). 
The radial equation of motion ([[]) and the continuity equation (Q) can then be written as a pair of 
ordinary differential equations for £ r and £h- 



dr 
dr 



- 2 - + A 

r 



dr J 
AP dlnP 



m?uj 2 rp n(n + 1) 



m 2 uj 2 p dr 



T X P 



fr 2 p 
TiP' 



Afr 

9 9 ' 



(10) 

(11) 



the solution of which requires two boundary conditions. At the surface of the star we take a free 
boundary: AP = 0, i.e. P' = —^ T dP/dr. The boundary condition at r = 0, where equations (|lO| ) 
and (11) have a regular singularity, is that the solutions be regular. Since at r ~ 0, P = P c + 0(r 2 ), 
p = p c + 0(r 2 ) and A = 0(r), where P c and p c are respectively the central pressure and density, 
this leads to £/, = £ r /n. 



2.1.2. Torque due to dissipation in the convective envelope 

The interaction between convective motions and the tidal flow is expected to lead to the 
dissipation of tidally excited waves (e.g. Zahn 1977). We model this effect as arising from a 
turbulent viscosity. To do this we suppose there is an additional dissipative force per unit mass 
acting in the convection zone given in spherical coordinates by: 



1 d ( 2 9v\ 



c pr 2 dr\dr) : ^ ^ 

where v is the flow velocity and v is the turbulent viscosity. Here, we assume that variation in the 
radial direction is the most significant and note that the viscous force is defined in such a way as 
to lead to a positive definite energy dissipation rate. 
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For the turbulent viscosity coefficient, we take (see, for example, Xiong, Cheng & Deng 1997) 



ci 



A 2 



t c 1 + c 2 (mt c /P ) 



(13) 



where ci and C2 are two constants, ci being on the order of unity, P a = 2tt/lo is the orbital period, 
A is the mixing length, and t c = l/\/\N' 2 \ is the convective timescale. The viscosity is then ciA 2 /t c 
for small forcing frequency muj. The factor 1 + (mt c / P ) s , where we shall use s = 2, allows for 
a reduction of efficiency of the damping of high frequency oscillations to which the convection 
cannot adjust. A similar prescription with s = 1 has been proposed by Zahn (1966) and used 
by Zahn (1989), whereas s = 2 has been considered by Goldreich & Keeley (1977) and used by 
Campbell & Papaloizou (1983) and Goldman & Mazeh (1991). Goodman & Oh (1997) have also 
recently put forward some arguments in favor of s = 2. For the mixing length, we shall take the 
standard relation A = a/ \dlnP/dr\ and set a = 3. 

In principle, equations (|l])-(||) should be solved with F c added on the right hand side of 
equation ([l]) in the convective envelope. However, this would increase the order of the differential 
system to be solved and make the numerical calculations much more complicated. Instead, we 
have found it adequate to solve first the adiabatic problem and then to treat the dissipative effect 
using a first-order perturbation theory. This is valid everywhere, except very close to resonances, 
because dissipation is weak enough so that the imaginary parts of £ r and ^ are much smaller in 
magnitude than their real parts. Thus we solve equations (ffl)-(S) without dissipative terms, and 
use these (real) solutions to calculate the rate of energy dissipation dE/dt due to convection which 
is given by: 



dE 
~dt 



[ pRe (F c ) • Re (v) dV, 



(14) 



where the integration is over the volume V c of the convective envelope and the angular dependence 
of F c and v has to be taken into account. Using the relation v = d£/dt and the expression (§) for 
£, we get 



dE 48?r 2 2 rRc 
dt 5 jr. 



pr 2 u 



dr ) V dr 



dr, 



(15) 



where Ri and R c are respectively the inner and outer radius of the convective envelope. Noting 
that the ratio of the rate of exchange of energy and angular momentum with the orbit is given by 
the pattern speed of the tidal disturbance, u, the torque exerted by the companion on the star is 
given by 



r 



1 dE 

lu dt 



(16) 
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This torque is positive because the star is non-rotating. 

When the frequency of the tidal wave is equal to that of some adiabatic normal mode 
frequency of the star, we can no longer use first-order perturbation theory because it would give 
an infinite torque. However, when the frequency u> is very close to a resonant frequency u>q, the 
torque will be given by an expression of the form 



~2f \2> 2' K J 

where A is an amplitude and 7 is the damping rate for the mode. First-order perturbation theory 
assumes dissipative effects are small in the response calculation. Therefore it is valid only for 
frequencies such that m 2 (uj — coq) 2 3> 7 2 . However, the damping rate, if small, can be found from 
first-order perturbation theory applied, as described above, very close to resonance where the 
mode dominates the response. Then it is given by (see, for example, Goldstein 1980) 



1 dE 

2 -< = ~m * (18) 



very close to resonance, where Ek is the kinetic energy of the mode: 



Ek = \J v P [Re(v]} 2 dV = 2 -^mW jf* pr 2 (g + 6$) dr, (19) 

the integration being over the volume V of the star (in the first integral, the angular dependence 
of v has to be taken into account). In fllg|), the total energy of the mode is 2Ek because at 
resonance there is equipartition between kinetic and potential energy. We calculate 7 using ( |l8|) as 
outlined above making sure that u is close enough to ujq by checking that 7 remains approximately 
constant when u is slightly changed. To get A, we fit the torque obtained from first-order 
perturbation theory to the expression ([17]) in the region approaching the resonance where still 
m 2 (io — ujq) 2 3> 7 2 . This procedure works satisfactorily when 7 is small with consequent strong 
resonances. This appears to be the situation when turbulent viscosity alone is assumed to act. 
However, radiative damping cannot be neglected for resonances at low forcing frequency and this 
is discussed below. 



2.1.3. Torque due to non-adiabaticity in the radiative core 

Non-adiabatic effects become important when the radiative diffusion time across the length 
scale associated with the tidal response shortens to become comparable to the wave propagation 
time across the star. In principle, these effects should be taken into account both in the radiative 
core and above the convection zone. However, since g modes, which are excited in the radiative 
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core, are evanescent in the convective envelope, we do not expect non-adiabaticity to play an 
important role above the convection zone. To take into account non-adiabatic effects, equation (|3|) 
has to be modified in the radiative core to: 

„T^) = _V.F<, (20) 

where F' is the perturbed radiative flux. The radiative flux is given by the radiative diffusion 
equation 

F = -KVT, 

where T is the temperature and K = AacT s / (Sup) is the radiative conductivity, with a being the 
Stefan-Boltzmann radiation constant, c the velocity of light and k the opacity. Therefore 



r 2 dr 



r 2 ( K— + K'— 

or dr 



V\ (KT>) , (21) 



where is the non-radial component of the operator V. We now suppose that close to 
resonance, the response behaves exactly like a free g mode with very large radial wavenumber k r 
so that WKB theory can be used together with the local dispersion relation to evaluate V • F'. 
The dominant term in the right-hand side of (|2l| ) is then —Kd 2 T'/dr 2 = Kk 2 T', and all the other 
terms can be neglected. For a high order free g mode, the local dispersion relation gives (see, for 
example, Unno et al. 1989): 

,2 _ iV 2 1(1 + 1) 

This expression of k r is derived under the adiabatic approximation. However, since we want to 
incorporate non-adiabatic effects to the lowest order, we do not need to take them into account 
in evaluating k r . We have k r 3> 1/r in the radiative core because N 3> mto there. The perturbed 
temperature T" can be expressed as a function of P' and p' using the thermodynamic relation 

¥-^-5-^. (^) 

T XT P XT P 
where xt = dlnP/dlnT) and xp = dlnP / dlnp) T . 

The equation of state in the radiative core of a solar-type star is primarly that of a perfect gas 
(we neglect here the radiation pressure which is very small compared to the gas pressure). We 
then have Ti = T3 = 5/3 and Xp = XT = 1- Using the above, the fact that I = 2, d/dt = —imiv 
and the relation (El), we can recast (pH) under the form 
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p l + ie TiP 1 + ie 
where 



1 + !£Fl P ' A fr, (24) 



16 - T4JV2 (25) 



5 (muj) KpPr 2 

For high order free g modes, we have (see, for example, Unno et al. 1989) 

P'/P moj dlnP 
Ai r ~ iV d/nr ' 

which means that P'/P <C A£ r in the radiative core of a solar-type star at the frequencies of 
interest. The non-adiabatic correction of the term associated with pressure in equation (p4) is 
then very small compared to the non-adiabatic correction of the term involving £ r . Therefore 
equation (p4| ) can be approximated by 

d P' 

- = — -A£ r , (26) 

p r 1 p ^ ' v ; 

where we have defined A = A/(l + ie). We note that because we have identified the tidal response 
with a normal mode, this calculation is valid only close to resonances. Equation ( p6[ ) is similar 
to (H) but with A being replaced by A. The system of differential equations we have to solve to 
get the non-adiabatic response of the star is then the same as (|lO|)-(|Tl|) but with the following 
modification in the radiative core. In equations ( |io|) and (11), A, where it appears as a coefficient 



of £ r , has to be replaced by A. The system of differential equations so obtained is complex. In 
general, for the periods we consider, e < 5 x 10~ 4 is small. We then calculate both the real and 
imaginary parts of the response, so that the torque can be computed directly from: 

T -I v "' a -W dv (27> 

where the integral is over the volume V of the star. The angular dependence of *$>t and p' has to 
be taken into account in this expression. Equation ( f27| ) can be recast under the form (Savonije & 
Papaloizou 1983): 

T = [ R& Im(p'ydr. (28) 



As in § gXg , we can calculate the damping rate 7' in resonances due to non-adiabatic effects 



using (|i~8[). Here, the rate of energy dissipation is calculated from the torque (see 16). 
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2.2. Low frequency limit 

In the limit of small \u>\, the following relations may be obtained for the adiabatic equilibrium 



tide: 



PL = ~fr 2 p, (29) 



and, if iV 2 ^ 0: 



Zr,eq = fr 2 p(^) , (30) 

^ = ^^ri {r% ^ (31) 

where the 'eg' subscript denotes the equilibrium value. Using these displacements, the torque may 
be calculated using equations (|l5|) and (lilf). 



However, we comment that equation (p0|), which states that the Lagrangian perturbation to the 
pressure is zero, can only be derived in the adiabatic low |w| limit if the Brunt- Vaisala frequency 
is not zero (in practice, one also requires that the forcing period be short compared to the thermal 
timescale of the star, but the latter is so long that it can be assumed to be infinite in this context). 
Equations (|3C|)-(|3l|) do not apply in a finite region where strictly N 2 = 0. In that case the fluid is 
locally barotropic, and the displacement can be written as the gradient of a potential: 



£ = V [*(r)y n ,, 
The continuity equation then gives for low frequencies 



(32) 



r 2 dr 



2 d< S>\ 
^ dr J 



n(n + 1)/9<E> 



-A 



PLP fr 2 p 2 



eq 



riP 



(33) 



Equation (33) gives a second order differential equation for <£(r). This applies inside the region 
where iV 2 = 0. It is possible, using the two available boundary conditions for (|33"1), to match ^ r ^ eq 
given by (|30|) at the boundaries of such a region, but not in general £h,eq given by (|3T|). This 
means that there will tend to be a discontinuity in the tangential displacement at the boundaries 
for low frequencies. 



When |iV 2 | is not zero but very small, in particular small compared to 



2, ,2 



m uj 



which corresponds 



physically to the convective timescale being much longer than the forcing period, the tidal response 
more closely matches that given by (33) than that given by equations (30)-([H]). This feature 
causes a very slow convergence towards the low frequency limiting solution {equilibrium tide) 
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found here, as well as near discontinuous behavior near the inner convective envelope boundary. 



This is borne out by our numerical results (see § 3.2.1 and § 



2.3. Timescales 

2.3.1. Orbital evolution and stellar spin up timescales 

The torque, T, gives the rate at which angular momentum is transferred from the orbit to the 
star. We can then calculate a tidal evolution (decay) timescale of the circular orbit: 

t rb = - , (34) 

where /i = M p Mq/ \M p + Mq) is the reduced mass. In principle, the variation of the torque 
with ijj has to be taken into account for the total decay time to be calculated. However, since the 
torque increases as the companion spirals in, t OT \, is mainly determined by the initial position of 
the companion, and a good estimate can be obtained using the above formula. 

This exchange of angular momentum also results in the spin up of the star, the timescale of 
which is given by t sp = Iu/T (Savonije & Papaloizou 1983), with / being the stellar moment of 
inertia. Out of resonance, angular momentum deposition initially occurs mainly in the convective 
envelope where the turbulent viscosity acts (Goldreich & Nicholson 1989). It is then of interest to 
calculate the spin up timescale for the convective envelope alone, which is t sP)C = I c u/T, where 
I c is the moment of inertia of the convection zone. We note however that on the long timescales 
associated with tidal evolution, angular momentum may be redistributed throughout the star. 



2.3.2. Circularization 

In practice, we find that the torque arising from a companion in a circular orbit varies 
with frequency approximately oc uA This result can be used to relate the orbital circularization 
timescale to the initial orbital decay timescale provided that the eccentricity is not too large. In 
practice, both these timescales can be significantly longer than the spin up timescale of the star 
due to its relatively small moment of inertia. We should then consider that the star is synchronized 
with the orbit. 

The ratio between the orbital decay timescale and the circularization timescale t c i rc is 
found to be about 6 for the calculated frequency dependence of the circular orbit torque (see, 
for example, the expressions given in Savonije & Papaloizou 1983, 1984). This appears to be 
independent of whether the star is assumed to be synchronously rotating or non-rotating in that 
the circularization rate calculated assuming no rotation, as we do here, gives a reasonable estimate 
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in the synchronous case also. In addition, to evaluate t c i rc for an equal mass system, we have to 
take into account the reciprocal torque exerted by the primary on the companion. To do this, we 
take t c i rc to be proportional to a factor which is 1/2 when M p = Mq and 1 when M p <C Mq. To 
a reasonable approximation, we then get: 

tarc = 6(1 + Mp/M )' (35) 



2.4. Circularization timescale as a calibration of turbulent viscosity 

One of the purposes of this paper is to calculate the tidally induced velocities on the star. In 
order to do this, the processes responsible for dissipating the disturbance should be included as 
accurately as possible. An important dissipation process is that associated with turbulent friction 
in the convection zone. As this process has been suggested as being responsible for circularizing 
the orbits of main sequence binaries of sufficiently short period (see Mathieu 1994), we investigate 
whether reasonable assumptions about the behavior of turbulent viscosity can lead to the observed 
circularization rates. 



2.4.1. Background 

Zahn & Bouchet (1989) have investigated the pre-main sequence evolution of late-type 
binaries in which the stars are fully convective. The main conclusion of their work was that 
circularization occurs at the very beginning of the Hayashi phase, with hardly any decrease of the 
eccentricity on the main sequence. The cutoff period they predict is between 7.3 and 8.5 days. 
According to them, observations show a cutoff period around 8 days independent on the age of 
the star and are then in agreement with their results. 

However, we note that for this agreement to be reached, some observations had to be 
discarded. Those by Mayor Sz Mermilliod (1984), which indicated that the cutoff period of 
late-type binaries was at most 5.7 days, and those by Mathieu Sz Mazeh (1988), which showed 
that the cutoff period in the 4 Gyr old cluster M 67 was more than 10 days. In a review article, 
Mathieu (1994) (see also Mathieu 1992 and Mathieu et al. 1992) has confirmed that "cutoff periods 
increase with age, consistent with active main sequence tidal circularization". The pre-main 
sequence cutoff period is very likely to be 4.3 d (an upper limit being 6.4 d), and cutoff periods 
for solar-mass binaries are 7.05 d in the Pleiades (0.1 Gyr), 8.5 d in the Hyades (0.8 Gyr), 
12.4 d in M 67 (4 Gyr) and 18.7 d in the halo field (16 Gyr). We therefore conclude that active 
circularization does take place for main sequence binaries and that it is less efficient than proposed 
by Zahn Sz Bouchet (1989) on the pre-main sequence. 

Recently, Claret Sz Cunha (1997) have applied the formalism developed by Zahn (1989) to 
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different stellar models. They have computed the parameters which enter into the expression for 
the circularization timescale, which is based on treatment of the equilibrium tide, for a wide grid 
of stellar models as a function of mass and time. Their conclusion is that turbulent dissipation is 
too low by a factor 100-200 during the main sequence to fit the observed cutoff periods. 



3. Numerical results 



The calculations presented in the following section are applied to the standard solar model 
described by Christensen-Dalsgaard et al. (1996). 



3.1. Equilibrium tide 

The torque associated with the equilibrium tide was calculated as indicated in § As 
mentioned there, this calculation is only expected to apply at very low frequencies. For periods 
between 4.23 and 12.4 d, and c\ = C2 = 1, the calculated torque can be interpolated by the 
following power law: 

T (g.cm 2 /s 2 ) = 1.200 x 10 55 ( ^ M \ T ) u 4 - 08 . (36) 
V y ' / \M P + M Q J v ; 



3.2. Dynamical tide 

For the calculation of the dynamical tide, we solve numerically the differential equations ( |i~0| ) 
and ( |Tl| ) using a shooting method to an intermediate fitting point (Press et al. 1986). To evaluate 
non-adiabatic effects in the radiative core close to resonances, we modify these equations in 
the way described in § 2.1.3| . We define the dimensionless quantity x = r/R c , where R c is the 



outer radius of the convective envelope. With this notation, the equations are integrated from 
Xi n = 10~ 6 to x ou t = 1.00071256. The radiative core extends from x = to x ~ 0.7. The results 
presented below have been obtained with the values c% = ci = 1 and a = 3 (a being the ratio 
of the mixing length A to the pressure scale height). We discuss the effect of changing these 
parameters. 



3.2.1. Tidal response and velocity at the surface of the star 

For illustration purposes, we plot the horizontal and radial displacements induced in the star 
at orbital periods of P Q = 4.23 d and 8.46 d away from resonance in the adiabatic approximation. 
The first period is that inferred for the system 51 Pegasi. The spatial distribution of the real 
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parts of mu^r and muj^h are shown in Figure |j. These represent typical values of the radial and 
horizontal velocities, the maximum values being three and six times larger respectively. Since 
these quantities depend on the perturbing mass through the ratio M p /(M p + Mq), they have been 
represented in units of this factor. 

As expected, the stellar response shows oscillations between turning points near the center 
and the inner radius of the convection zone, where m 2 uj 2 = N 2 , otherwise it is evanescent. The 
horizontal displacement varies rapidly in the photosphere because the temperature drops to zero 
rapidly there. 

We see from Figure [l] that, when the perturbing mass M p = Mq, the maximum radial 
velocity at the surface of the star is about 6 and 1 m/s for P Q = 4.23 d and 8.46 d respectively. 
These numbers drop to 10 -2 and 2 x 10 -3 respectively when the perturbing mass is one Jupiter 
mass (M p = 1CT 3 Mq). 

The radial displacement and the perturbed pressure at the surface of the star are well 



approximated by the equilibrium values (30) and (29) respectively. These quantities are found 
to be insensitive to the existence of resonances with the consequence that the radial velocity at 
the surface of the star never differs much from the values given above. For the smallest periods 
considered, the ratio \Re{£,h) / Re{£, r )\ at the surface of the star can vary by up to one order of 
magnitude on passage through resonance. This is due to the fact that £/, is proportional to 
(P' - P^j jP (see |). Even though this ratio remains small, it can vary by up to an order of 
magnitude as a resonance is passed through. 

The numerical results indicate that both the amplitude and the wavelength of the response 
increase with the orbital frequency, in agreement with the theoretical expectation of a smaller 
radial order for higher frequencies (Christensen-Dalsgaard & Berthomieu 1991 and references 
therein). 



Finally, as expected (see § 22), the plots shown in Figure [T] (see also Figure |3|) indicate that 



at the boundary of the radiative core and the convection zone there is a near discontinuity in the 
mean value of £/j obtained after averaging out the interior oscillations. 



3.2.2. Circular orbit torque 

i) Resonances: 

Figure ^ shows T versus u in a log-log representation for three different small frequency 
intervals in the vicinity of P Q = 4.23, 8.5 and 12.4 d. On each plot, the dotted line gives the values 
obtained from the theory of the equilibrium tide as given by (^). Since the torque depends on 
the perturbing mass through the factor M 2 /(M p + Mq) 2 , it has been represented in units of this. 
These plots show several resonances, which occur when the frequency of the tidal wave is equal 
to the frequency of some normal mode of the star. The left panels show the torque arising from 
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convective dissipation, through turbulent viscosity, alone (see § 2.1.2). These plots have been 
displayed for comparison with the right panels, for which radiative damping has been taken into 
account in the resonances (see § |2.1.3j ) . 

As indicated by the Table [l], the normal mode damping rate due to radiative damping (7') is 
much larger than that due to convective dissipation (7). 



Table 1: Normal mode damping rates 



Po 




7 




7' 


(days) 


is' 1 ) 








4.23 


1.72 x 10" 5 


6 x 10" 


12 


1Q -io 


8.5 


8.56 x 1(T 6 


5 x 10" 


12 


7 x 10- 10 


12.4 


5.86 x Id -6 


4 x 10" 


12 


2 x 10" 9 



Thus the torque in the center of resonances, where they are significant, is predominantly 
determined by radiative damping. For the frequencies we consider, this contribution to the 
torque becomes much smaller than that due to turbulent viscosity in the tail of the resonances. 
This means that non-adiabatic effects in the radiative core are negligible away from resonances. 
Therefore, in the right panels of Figure [2|, we have plotted the torque resulting from radiative 
damping acting alone in the center of resonances, and that resulting from convective dissipation 
acting alone away from resonances. A comparison between the strength of the resonances shown 
in the right and left panels indicates the importance of non-adiabatic effects in the radiative core. 
As expected, the resonances are weakened and broadened, this effect being marginally important 
for P a ~ 4.23 d. 

We now discuss the properties and the effect of resonances on the tidal torque. From now on, 
when resonances are discussed, we shall refer to the calculations which take into account radiative 
damping. 

In the neighborhood of P Q = 4.23, 8.5 and 12.4 d, the relative separation, Aw/w, between 
successive resonances is respectively 4.5 x 10 -3 , 2 x 10 -3 and 10~ 3 . The relative width, Su/u, of 
the resonances is respectively 3 x 1CP 4 , 1.5 x 10 -4 and 10 -4 . Here we have arbitrarily defined the 
width of a resonance as being the frequency interval over which T is at least 3 times larger than 
the minimum torque obtained just out of this resonance. 

To calculate the probability of the companion being in a resonance, we have to take into account 
the fact that it drifts away from the resonances much more rapidly than elsewhere. The relevant 



quantity for calculating the tidal evolution timescale is 1/T (see expression 34). For a fixed 
oscillation spectrum, we can approximate this probability by 5uj/Aoj times the ratio of the mean 
value of 1/T over a resonance to the mean value of 1/T between two resonances, where the mean 
value is defined by 
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l \ _ r 
r/~ J ~r 



du , 



with the integrals being taken over the relevant frequency interval. 

This gives a probability of being in a resonance which is close to 0.7% for P Q = 4.23 and 8.5 d and 
2% for P Q = 12.4 d. The fact that the probability of being in a resonance increases with P Q is not 
significant, because resonances get weaker when the period increases (see Figure §). 

We note that this discussion applies only if the a priori probability of being in any frequency 
interval of a given width is independent of the frequency as might be expected to be a reasonable 
assumption if the normal mode spectrum is fixed. However, different circumstances may apply 
if the combined effect of orbital and stellar evolution were to lock the companion in a resonance 
with changing location. But we shall not consider the possibility of this process here. 

As the companion spirals inwards, it goes through a succession of resonances. However, for a fixed 
normal mode spectrum, the above calculation tells us that its migration is controlled essentially 
by the non-resonant interaction. This can be seen by comparing < 1/T > evaluated over a large 
frequency range, both taking into account and neglecting the resonances. Such a comparison 
shows that neglecting resonances changes < 1/T > by at most a few per cent. 

ii) Relation between the mean torque, to and the circularization timescale: 

Here we interpolate the numerical results to express the torque as a power of the frequency. 
To begin with, we consider the three frequency intervals described above. We take the appropriate 
torque to be 1/ < 1/T > where the mean values are taken over the frequency intervals displayed 
in Figure H The results can be interpolated with the following relation: 

T (g.cm 2 /s 2 ) = 1.654 x 10 53 ( — i%— V c 3 - 85 . (37) 
^ 1 > \M p + Mq) v ; 

We have checked that the above formula gives a good estimate of the torque at other non-resonant 
frequencies between 4.23 and 12.4 d. 



Since the index of the power law ( |37| ) is close to 4, the circularization timescale t c i rc is given 
by (H). At P = 12 A d, t c i rc is found from the above formula to be 56 times larger than the 



timescale of 4 Gyr that is indicated by the observations. 

We note that both the dynamical and equilibrium tide calculations give a power law with an index 
close to 4, which results in the circularization timescale being proportional to the binary period 
P Q to the 13/3. For comparison, Zahn (1977), Zahn (1989) and Goldman & Mazeh (1991), using 
equilibrium tide calculations, found t c i rc to be proportional to P a raised to the power 16/3, 13/3 
and 10/3 respectively. The difference between these results can be related to a different choice of 
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s in the expression (|T^) for v. These authors used an expression similar to ( |i~3| ) with respectively 
s = 0, 1 and 2. The fact that we obtain an index close to 13/3 by setting s = 2 or even s = 1 (see 
below), in contrast to the results above, is at least partially due to the effectively smaller value of 
C2 we used (see below). 

We comment further that Tassoul (1988) found t c i rc oc Po 9 ^ 12 for his postulated alternative 
hydrodynamical mechanism for tidal circularization. 



3.3. Comparison between calculations based on the dynamical and equilibrium tides 

The results presented above show that the torque corresponding to the dynamical tide is 
smaller than that given by the adiabatic equilibrium tide for all the frequencies we have computed. 
However, the difference tends to decrease as the frequency gets smaller. From expressions (j36|) 



and (37), we calculate that the ratio of these torques is indeed about 6.0 and 4.8 for P Q = 4.23 and 
12.4 d respectively. 

Figure || shows ^ r>eq /R c , £, r /Rc and £,h/Rc in units M p /(M p + Mq) versus x in the range 
0.6 < x < x out for P Q = 4.23 and 12.4 d. 

It is clear from these plots that £ r departs from the asymptotic value in the convective 
envelope. The difference is not large, being about 17%, but the derivatives of £ r and from which 
the torque is calculated (see |l~5|) depart more from their asymptotic values. 



In the limit where the magnitude of the Brunt- Vaisala frequency is everywhere large compared 
to the tidal forcing frequency, calculations based on the dynamical tide should converge towards 
those based on the equilibrium tide. This is because the convective timescale is small enough that 
the convective motions adjust essentially instantaneously to the tidal forcing. 

We have checked this expectation by artificially increasing the magnitude of the Brunt- Vaisala 
frequency in the convection zone. Except in the part of the convective envelope just below its outer 
radius, wherever |iV 2 | < quj 2 , q being an arbitrary constant, we make the replacement N 2 = —qui 2 . 
In Figure |4] we plot £, r ,eq/£,r versus x in the range 0.6 < x < x ou t for P a = 12.4 d and for q = 10, 
100 and 400. For comparison we also plot the case corresponding to the original solar model. 

As expected, C r ,eq/Cr converges towards 1 when the magnitude of the Brunt- Vaisala frequency 
is increased. We note that q = 400, which corresponds to \N\ > lOmuj, gives 1 — £r,ecj/£r 
smaller than 5%. This confirms that the asymptotic limit is reached when the magnitude of the 
Brunt- Vaisala frequency is very large compared to the frequency of the tidal wave. 

We note that the very slow convergence towards the equilibrium tide was predicted from 



the arguments presented in § 2.2. Also expected was the discontinuity in the mean value of 
obtained after averaging out the interior oscillations, at the boundary of the radiative core and the 
convection zone that is observed in Figure ||| (see also Figure [l]) . 
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3.4. Calibration of the turbulent viscosity 

We shall limit the comparison of our results with observations of main sequence binaries 
because our calculations do not apply to pre-main sequence stars, which have a much larger 
convective envelope than the Sun. As mentioned above, the observed circularization timescale 
we have to fit is then 4 Gyr for P Q = 12 A d (Mathieu 1994). As indicated above, when using a 
simple estimate of the turbulent viscosity based on mixing length theory for non-rotating stars, 
the circularization timescale we get from our calculations for this period is 56 times larger than 
4 Gyr. This indicates that either 

(i) solar-type binaries are not circularized through turbulent viscosity acting on tidal perturbations 
(but see Tassoul 1988 and Kumar & Goodman 1996 for other suggested tidal mechanisms), or 

(ii) dissipation in the convective envelope of solar-like stars is significantly more efficient than is 
currently estimated. 

Tassoul (1995) postulates that efficient tidal dissipation occurs in a very thin Eckman layer 
close to the surface of a tidally deformed star and that this process greatly increases the efficiency 
of tidal interactions. But a refutation of the notion that the free surface boundary condition 
appropriate to the tidally deformed star, rather than the more common rigid boundary condition, 
leads to such an effective boundary layer, has been given by Rieutord & Zahn (1997). Further 
Tassoul & Tassoul (1997) state that their mechanism is inapplicable to extreme mass ratio cases 
such as 51 Pegasi that we consider later in the paper. 

We now consider briefly here the mechanism proposed by Kumar & Goodman (1996), namely 
enhanced dissipation associated with high order oscillation modes excited through parametric 
instability. The growth rate for the most rapidly growing modes is expected to be a ~ mu;^ r / Rq, 
where £ r is the radial displacement in the primary oscillation, which we shall assume to be the 
equilibrium tide, evaluated at r = Rq, and muj is its frequency. 

If we assume that the non-linear development of the parametric instability and subsequent 
dissipation of the excited modes leads to an effective viscosity, and frictional dissipation rate 
tj 1 = a, which is big enough to suppress the linear instability, then we expect from the classical 
theory of Darwin (1879) that there will be a phase lag 9t associated with the tide given by 



For a binary of unit mass ratio and period ~ 10 d, synchronization occurs on a timescale 
very much shorter than that required for circularization, so that we assume that the stellar 
rotation is synchronized with the orbit and m = 1 in the calculation of a. For small eccentricity, 
the circularization timescale t c i rc is approximately given by l/t c i rc = £l a 9t, where the apsidal 
motion frequency is Q a = 15k(M p / 'Mq)(Rq / D) 5 u , with k being the apsidal motion constant 
(Cowling 1938). We use the equilibrium tide value ( |30|) to estimate ^ t /Rq at r = Rq as 
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(Rq/D) 3 /4 for a mass ratio of unity. We then get 1/Wc = 7 '.bkuj{RQ / 'D) u oc w 25 / 3 for 9 t given 
by (|38|). For and orbital period of 12.4 d, this gives t c i rc = 6 x 10 3 //c Gyr. Since k ~ 0.01, it is not 
very likely that this mechanism will be able to explain the observed circularization rates. However, 
we stress that this has not been shown from a full non-linear calculation of the development of 
parametric instability. 



3.4-1- Enhanced turbulent viscosity 

In general, a large increase in the simply estimated turbulent viscosity coefficient is needed in 
order to explain the observed circularization rate. We now investigate what is needed to achieve 
this and give the numerical results of tests we have carried out. In all cases, we have checked that 
the velocity at the surface of the star is not sensitive to the magnitude of the turbulent viscosity 
assumed. 

In addition, as indicated above, the resonances are essentially controlled by radiative damping 
in the radiative core as long as 7' S> 7- When the turbulent viscosity is enhanced, 7 is increased. 
However, in the tests we present below, for orbital periods larger than ~ 8 days, 7' stays large 
enough compared to 7 so that, although the non resonant torques increase, the central structure 
and strength of the resonances is determined by radiative damping. For orbital periods on the 



order of ~ 4 days, 7 can become comparable to 7'. Then the damping factor in 17 has to be 
replaced by (7 + j') 2 , and the strength of the resonance is reduced by a factor ~ 4. 

We first consider the effect on the circularization timescale of varying the parameters c%, c 2 , 



s and A in the expression (13) for v. Note that the denominator we used in this expression 
is 1 + ci(mt c j P ) 2 rather than 1 + c 2 (mu;i c ) 2 . Using the latter with C2 = 1 is equivalent to 
setting c 2 = (27r) 2 in the former. At present, it seems that our knowledge of convection does not 
allow discrimination between these possibilities (see, for example, the discussion in Zahn 1989). 
However, we note that Goodman & Oh (1997) have recently put forward some arguments in favor 
of c 2 = (2vr) 2 . Using c 2 = (2vr) 2 in @ results in a circularization time for P a = 12.4 d that is 10 
times larger than that obtained with c 2 = 1. If we set c 2 = 0, t c i rc is decreased by only a factor 
2. This is because when c 2 = 1 the factor c 2 (mt c /P ) 2 is already smaller than, or even very small 
compared to, unity in a large part of the convective envelope. Thus it seems that adjusting the 
way in which turbulent viscosity responds to short period forcing cannot produce the required 
enhancement in this case. 

In Table || we summarize the results obtained for different values of s and c 2 , and in Table ^ 
we indicate the corresponding index of the power law in expression (|3~7"|). We note that c 2 = 1 
corresponds to the denominator in ( |i~3|) being 1 + (mt c / P Q ) S , whereas c 2 = (2ir) s corresponds to 
1 + (mujt c ) s . 



We note that setting c 2 = (2-7r) 2 with s = 2 decreases the index of the power law in (|37|) down to 
~ 3.3. This gives t c i rc proportional to the orbital period to the ~ 11/3, which is similar to the 
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Table 2: t c j rc for different values of s and C2 



Po 


C2 


s 


tcirc 


(days) 






(Gyr) 


4.23 


1 


2 


2.46 






1 


1.37 


12.4 




2 


220 






1 


211 


4.23 


(27T) 2 


2 


39.4 




2tt 


1 


5.98 


12.4 


(27T) 2 


2 


2010 




2tt 


1 


685 



Table 3: Index of the power law in (|37| ) for different values of s and C2 



C2 s Index of the power law 



1 


2 


3.85 




1 


4.3 




2 


3.3 


2tt 


1 


4.1 



value found by Goldman & Mazeh (1991). 

Zahn & Bouchet (1989) have argued that the prescription s = 2 suggested by Goldreich & 
Keeley (1977) and used later by Campbell & Papaloizou (1983) and Goldman & Mazeh (1991) 
(see also Goodman & Oh 1997) would lead to cutoff periods in clear conflict with the observational 
data which they claim require s = 1. The results presented for the model adopted here do not 
support this statement. If C2 = 1, t c i rc hardly changes when s is changed from 2 to 1. If C2 = (27r) s , 
the difference between s = 1 and s = 2 is not dramatic for P Q ~ 8 — 12 d. Within the uncertainties 
associated with convection, our results do not allow a distinction to be made between s = 1 and 
s = 2. 

Finally, taking the mixing length to be the distance to the top boundary of the convective envelope 
rather than 3 times the pressure scale height does not affect significantly the circularization 
timescale. 

We note that the torque is directly proportional to c\ and a 2 . However, c\ is expected to be on 
the order of unity, and a is usually taken to be between 1 and 4. 
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3.4-2. Modifications to the Brunt-Vdisdld frequency 

The magnitude of the turbulent viscosity given by ( |i~3"D would be increased if the magnitude 
of the Brunt- Vaisala frequency was larger in the convection zone. This is because the convective 
timescale decreases. We comment that a reduction in the convective timescale, while maintaining 
the same length scale, implies larger convective velocities which would have to occur without 
increasing the heat flux. This is the essential feature of the modification. To illustrate the effect 
of increasing the Brunt- Vaisala frequency we consider the following dimensionless number: 

_ N 2 fdlnP 
77 9 V dr 

which is the superadiabatic temperature gradient when radiation pressure and variations of the 
mean molecular weight are neglected. The accuracy with which rj is known from helioseismic 
observations is not better than ~ 10~ 2 (Gough 1984). However, in most of the convective envelope, 
this parameter, estimated from mixing length theory applied to a non-rotating stars, is much 
smaller than 10 -2 . 

We have made a numerical investigation in which we increased |iV 2 | in the convection zone 
by replacing rj by min(p?7, 10~ 3 ), p being an arbitrary constant, wherever rj < 10~ 3 (except just 
below the outer radius of the convective envelope). We have considered p = 50 and p = 100. It is 
doubtful that consideration of present helioseismic data could preclude such an increase of | iV 2 1 
(Thompson 1997). 

Figure [| shows rj in the convective envelope versus x for the original solar model and for p = 50 
and p = 100. We also display the factor by which |iV 2 | has been increased in each case. 

The circularization timescales we find with C2 = 1 and s = 2 for P Q = 12.4 d when p = 50 and 
p = 100 are respectively 13 and 7.7 Gyr, which are now larger than the observed value by factors 
of 3 and 2 respectively. Small discrepancies of this magnitude could be dealt with by adjustments 
to the mixing length or c\. When either p = 50 or p = 100, the circular orbit torque is found to be 
proportional to a; 4 ' 6 . 

If, keeping s = 2, we adopt C2 = (27r) 2 , t c i rc is increased only by a factor 1.3 compared to the 
case C2 = 1 for p = 100. This is because the factor C2(mt c /P ) 2 is smaller than, or even very small 
compared to, unity in almost all the convection zone whatever C2 between 1 and (27r) 2 . 

Although an increase in |iV 2 | of the magnitude we consider might be thought to be unrealistic, 
we note that such an increase in the deep layers of the convection zone has also been considered 
by D'Silva (1995) as a means of explaining the dynamics of sunspots without invoking too strong 
a magnetic field. In his model, which applies strictly to a star rotating at the same rate as the 
sun, \N 2 \ has to be larger than 4 x 10 -11 s -2 , which means that it has to be multiplied on average 
by a factor ~ 8. In the solar model we use, we need to multiply |iV 2 | by a factor between 100 
and 400 for 0.722 > x > 0.713, between 10 and 100 for 0.828 > x > 0.722 and between 1 and 10 
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for 0.915 > x > 0.828 in order to get such a minimum value. If we do this, the circularization 
timescale we get for P Q = 12.4 d is only 6 times larger than the observed one. In this context, it is 
possible that if the proposed increase in |iV 2 | is related to the stellar rotation, this may be even 
greater for the more rapidly synchronously rotating star that is expected in the equal mass binary 
case. 

We note that numerical simulations of turbulent convection in the presence of rotation show 
an increase of |iV 2 | with the effect of rotation (Brummell, Hurlburt & Toomre 1996). This is 
because, as pointed out by Brummel et al. (1996), rotation influences the thermodynamic mixing 
properties of the convection in such a way that it leads to a decrease in correlation between 
temperature fluctuations and vertical velocities. The efficiency of the vertical convective transport 
is then weakened, with a subsequent enhanced superadiabatic mean stratification in the interior 
of the fluid (see their Figure 8. a). 

This suggests that the magnitude of the Brunt- Vaisala frequency in the convective envelope of 
rotating stars is actually larger than the values given by the solar model we have been using 
here. However, it seems questionable that the extremely large increase required to account for the 
observed circularization rates can be achieved. 



3.4-3. Turbulent Viscosity below the convection zone 

Another means of increasing the total amount of dissipation is to assume that turbulent 
viscosity acts down to some depth below the inner boundary of the convective envelope. This 
might be expected if convective overshooting takes place. However, this might not be a very 
effective process because of the slow convective motions expected and the rapid increase in |iV 2 | 
that occurs as the radiative zone is entered. 

Estimates based on the observed solar oscillation frequencies give an upper limit between 0.05 
(Basu 1997) and ~ 0.1 (Christensen-Dalsgaard, Monteiro & Thompson 1995) times the pressure 
scale height on the extent of overshoot below the convection zone. 

Here we consider a simple illustrative situation in which convective blobs or some other turbulent 
motions are able to penetrate into the stratified radiative core over some fraction z of the pressure 
scale height \dlnPfdr\~ 1 producing a turbulent viscosity. We model this by setting v to be 
constant from a distance 0.5z \dlnPfdr\~ 1 above the inner radius of the convective envelope down 
to the same distance below this radius, equal to its value at the top of this zone. Because such a 
viscosity is able to act on the short wavelength part of the tidal response associated with g modes 
it has a dramatic effect. 

For z = 1 and c\ = c<i = 1, the circularization timescale for P Q = 12.4 d is decreased by a factor 
about 400, being now 8 times smaller than the observed timescale. If we set C2 = (27r) 2 (see 
discussion above), i c j rc is increased by a factor 15, being about 4 times larger than the observed 
timescale. 
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For z = 0.1 and c\ = C2 = 1, we get a circularization timescale for P Q = 12.4 d about 30 
times larger than the observed one. For the calculated timescale to be in agreement with the 
observations, we need z between 0.4 and 0.5 with c\ = C2 = 1. In the model we have adopted, the 
effect is of less importance for shorter periods because the number of oscillations of the response 
in the region of the radiative core where turbulent dissipation is introduced decreases with forcing 
frequency. Therefore, the torque does not vary with w as a simple power law with an index close 
to 4. However, in reality, the effectiveness of the turbulent viscosity should be reduced for short 
wavelength disturbances giving a compensating effect to make it relatively less effective at low 
frequencies. 

The above calculations show that overshooting is not likely to be efficient enough to decrease 
the circularization timescales by a factor of about 50. To get such an effect, we indeed require an 
extent of overshoot below the convection zone at least 5 or 10 times larger than that deduced from 
the observations. In addition, we have not taken into account the fact that overshooting leads to 
an increase of the g modes length scale through a decrease in the buoyancy or the magnitude of 
N 2 . This in turn would decrease the amount of turbulent dissipation associated with overshooting. 



3.5. Fitting the observations 

As we have already mentioned above, calculations based on both the dynamical and equilibrium 
tide theories give a torque proportional to the orbital frequency raised to a power ~ 4 (see [36] 
and 37). If circularization of solar-type binaries does occur through the action of turbulent 



viscosity on the tides, then its magnitude has to be calibrated so as to account for the observed 
timescale. We have discussed in the previous section, somewhat speculatively, how the required 
enhancement of the magnitude of the viscosity above that obtained from simple estimates might 
be envisaged to occur. Since the enhancement might depend on forcing frequency, it is not clear 
that the resulting torque will still be proportional to the frequency to the ~ 4. However, the 
increase to |iV 2 | described above gave torques that approximately preserved this power law so that 
in the absence of additional information, we shall suppose it holds. Then the calibration acts only 
to adjust the coefficient of the power law. 

We note that the observations do not rule out any exponent between 3 and 5 (see below). Since 
our calculations can only strictly be applied to solar-type stars, we calibrate our results using 
tdrc = 4 Gyr for P Q = 12 A d. This gives (in cgs): 

T (g.cm 2 /s 2 ) = 5.086 x 10 35 ( ) ( ^-rV, (39) 

v y ' ) \m p + MqI yw^s- 1 ) 



or, equivalently: 
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r( 9 . Cm >^ 1.423 x!0-( ] ^_) 2 (A)- 4 . {40) 
The corresponding formulae for the orbital and spin up timescales (in giga-years) are: 

(m p /m + 1) 5/3 /f^y 3 ' 3 



^' G ^ = 2 - 763 xl ° mJm b [u) • < 41 > 

and t<j PjC = I c t sp /I (for the solar model we use, / = 1.064 x 10 54 g.cm 2 and I c = 1.5 x 10 53 g.cm?). 

Since the torque is proportional to u 4 , we can use the relation (^) for t c i rc: so that the 
circularization time is given by 

Even though our calculations can only be applied to solar-type stars, it is of interest to compare 
the circularization timescales we get from ( |4"3| ) with the observed ones. For P Q = 4.3, 7.05, 8.5 and 
18.7 d, (H) gives respectively t c i rc = 0.04, 0.3, 0.8 and 24 Gyr, to be compared with the observed 
timescales 0.003, 0.1, 0.8 and 16 Gyr respectively. The agreement for P Q > 8.5 d is within a factor 
1.5. For smaller periods, circularization is expected to occur when the convective envelopes of the 
stars are larger, making turbulent dissipation more efficient. We note that Mathieu et al. (1992) 

13 /3 

have already pointed out that a power law t c i rc oc P provides a close fit to the slope of the 
observed cutoff periods. However, the observations are equally well fitted with an index of 10/3 
and an index of 16/3 cannot be ruled out (Mathieu et al. 1992). 

If the simple estimate of the turbulent viscosity based on mixing length theory for non-rotating 
stars is used, the coefficients in the formulas for the torque have to be divided by a factor ~ 50 
whereas those in the formula? for the timescales have to be multiplied by the same factor. 

In Figure ||, we have plotted t or h in units (M P /M Q + 1) 5/3 / (M p /Mq) and t circ in units 
(M p /Mq + l) 2//3 / (M p /Mq) versus to and P Q in a log-log representation. 



4. Discussion 
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4.1. Application to 51 Pegasi 

It is of interest to apply these results to the system 51 Pegasi, for which the orbital period 
(assuming the observed oscillations are due to a companion) is P = 4.23 d. If the companion is 
a Jupiter mass planet (M p = 10 _3 Mq), then the tidal orbital evolution timescale given by (fll|) 
is t or b ~ 140 Gyr, the star spin up timescale (f42[) is t sp ~ 130 Gyr and the spin up timescale of 
the convective envelope is t sp ^ c ~ 18 Gyr. All of these timescales are long compared with the 
inferred age of 51 Pegasi (Edvardsson et al. 1993). If the companion is a low-mass star of O.IMq, 
as has been recently suggested, t or i, is 100 times smaller while t sp and t sPtC are 10 times smaller. 
We then expect the primary star to be synchronized with the orbit, in which case exchange of 
angular momentum is no longer taking place. Synchronization is actually expected if the mass of 
the companion is larger than about 10 Jupiter masses. The orbital decay timescale is also smaller 
than the age of the system, but since t sp < i or &, tidal interaction stops before the companion has 
plunged into the central star. 

If the simple estimate of the turbulent viscosity based on mixing length theory for non-rotating 
stars is used, all these timescales have to multiplied by ~ 50. In that case, synchronization is 
expected if the mass of the companion is larger than about 70 Jupiter masses. 

The planetary companion interpretation has been questioned recently by the reported 4.23 d 
modulation in the line profile of 51 Pegasi (Gray 1997), and the possibility that this modulation 
may be due to g mode oscillations has been considered (Gray & Hatzes 1997). 

We note that, according to our results, such a modulation could not be due to g mode 
oscillations tidally driven by a companion. For the oscillation to have a period of 4.23 d, the 
orbital period would have to be 8.46 d. The maximum perturbed radial velocity at the surface of 
the star induced by the companion would then be between 2 x 10~ 3 and 1 m/s for a perturbing 
mass between 10 -3 and 1 solar mass. These numbers do not depend on the magnitude of the 
turbulent viscosity assumed, and are not expected to be affected by the possibility of resonance. 
These velocities are at least about 50 times smaller than the observed ones. 

4.2. Summary 

In this paper, we have studied the dynamical response of the star to the tidal perturbation of 
a companion. We have computed the torque due to dissipation in the convective envelope using 
first-order perturbation theory. In the vicinity of resonance, we have also calculated the torque 
due to non-adiabaticity in the radiative core using a WKB treatment. We have found that the 
torque at effective resonances is mainly determined by radiative damping. We have carried out an 
analysis based on the adiabatic equilibrium tide, and showed that agreement with the dynamical 
tide calculations can be rather poor. For the unmodified stellar model and the periods of interest 
of several days, the torque derived using the equilibrium tide is 4 to 6 times larger than that 
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corresponding to the dynamical tide. 

We have found that the presence of fixed resonances do not affect the long term orbital 
evolution of the binary, so that the different timescales (orbital evolution, circularization and 
spin up) are mainly determined by the non-resonant interaction. Our calculations show that the 
viscosity that is required to provide the observed circularization rates of solar-type binaries is 
~ 50 times larger than that simply estimated from mixing length theory for non-rotating stars. 

We have explored some means by which this viscosity might be enhanced. We have found 
that it could become large enough if the magnitude of the Brunt- Vaisala frequency in the deep 
convective envelope were increased sufficiently. Such an increase is expected to be produced by the 
effect of rotation on convection, but it is questionable whether it can be of sufficient magnitude. 

We note that the strength of the resonances for orbital periods larger than ~ 8 days and 
the perturbed velocity at the surface of the star are insensitive to the magnitude of the turbulent 
viscosity assumed. Only for periods ~ 4 days is the strength of the resonances decreased by a 
factor ~ 4. The effective widths of the resonances affecting the tidal torques as well are reduced 
when the viscosity is increased. 

We have applied our results to 51 Pegasi, and showed that the oscillations which have been 
observed at the surface of this star cannot be a tidally driven non-radial g mode. Also we 
have found that the stellar rotation and the orbital motion of this system are expected to be 
synchronized if the mass of the companion exceeds 0.1 solar mass. 

We are grateful to M.J. Thompson for supplying us a solar model, and thank D.O. 
Gough, R.D. Mathieu and M.J. Thompson for helpful discussions. This work is supported by 
PPARC through grant GR/H/09454 and by NSF and NASA through grants AST 93-15578 and 
NAG 5-4277. C.T. and D.N.C.L. acknowledge support by the Center for Star Formation Studies 
at NASA/ Ames Research Center and the University of California at Berkeley and Santa-Cruz. 
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Fig. 1. — Real part of muj^ r (solid lines) and mw^ (dotted lines) in units M p /(M p + Mq) m/s 
versus x for xi n < x < 0.01 (top panels), 0.01 < x < 0.1 (middle panels) and 0.1 < x < x ou t 
(bottom panels), and for P Q = 4.23 d (left panels) and P Q = 8.46 d (right panels). These represent 
typical values of the radial and horizontal velocities, the maximum values being three and six times 
larger respectively. 
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Fig. 2.— log 10 {T) with T in units M 2 /(M p + M ) 2 g.cm 2 /s 2 versus %i (w) for P D = 4.23 d 
(top panel), 8.5 ci (middle panel) and 12.4 d (bottom panel). The solid and dotted lines correspond 
respectively to the dynamical and equilibrium tides calculations. On the left panels, T is calculated 
using convective dissipation only. On the right panels, T is calculated using radiative dissipation 
alone in the resonances and convective dissipation alone away from resonances. 
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Fig. 3. — Cr,eq/Rc (solid line), £ r /Rc (dashed line) and £h/Rc (dotted line) in units M p /(M P + Mq) 
versus x for 0.6 < x < x out , and for P Q = 4.23 d (top panel) and 12.4 d (bottom panel). 
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Fig. 4. — Cr,eq/Cr versus x in the range 0.6 < x < x ou t for P Q = 12 A d and for q = 10 (dotted line), 
100 (short-dashed line) and 400 (long-dashed line), the definition of which is given in the text. For 
comparison we have also plotted the case corresponding to the original solar model (solid line). 
(,r,eq/£,r is close to 1 when |iV| muj. 




Fig. 5. — Top panel: Dimensionless parameter n = N 2 /(g dlnP/dr) in the convective envelope 
versus x. The curves correspond to the original solar model (solid line), and to the models with 
increased |iV 2 | (p = 50, dotted line, and p = 100, dashed line, where p is defined in the text). 
Bottom panel: Factor by wich |iV 2 | is increased in the convective envelope versus x. inal 
corresponds to the original solar model, and iV 2 corresponds to the models with p = 50 (dotted 
line) and p = 100 (dashed line). 
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Fig. 6.— t orb in units (M p /Mq + l) 5 ' 3 / (M p /Mq) Gyr (expression gTJ) and t circ in units 
(M p /M Q + 1) 2/3 / (Mp/Mq) Gyr (expression H) versus w and P in a log-log representation. The 
cross indicates the position of 51 Pegasi. These timescales fit the observations. // instead they are 
calculated using the simple estimate of the turbulent viscosity based on mixing length theory, they 
have to be multiplied by ~ 50. 



